Tunneling anomaly of superconducting films in strong magnetic fields. 
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We consider the tunneling Density of States (DoS) of superconducting films driven to the param- 
agnetic phase by the Zeeman splitting. We show that there is minimum in the DoS whose position 
depends on the orientation of the applied field. This dependence, not predicted by previous theo- 
retical calculations, is in agreement with a recent experiment. 



It is well known that superconductivity can be de- 
stroyed by applying a magnetic field because of the break- 
ing of the time reversal symmetry (see e.g. Ref. 0). The 
magnetic field acts on both the orbital motion and the 
spin of the electrons; while in the bulk the orbital effect 
is responsible for the suppression of superconductivity, in 
low-dimensionality systems (d < 2) the Zeeman splitting 
of the states with opposite spin and which are degen- 
erate in zero field can be the dominant mechanism for 
this suppression^^ In the normal phase, the theory of 
interaction corrections to the Density of States (DoS) v 
reviewed in Ref. predicts the appearance of singular 
contributions to the DoS located at the Zeeman energy 
Ez = gL/J-BH (gL is the Lande g-factor, fiB the Bohr 
magneton and H the magnetic field). However experi- 
ments performed almost a decade ago£ and subsequent 
theoretical calculations 3 showed that there are more sin- 
gular corrections located at a lower energy E + : 



E+ = {E z + fl) /2 , 



where 



n = \/e 2 z - a 2 
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and A is the BCS gap. In a small grain, the origin of these 
corrections can be understood as follows^ when a spin- 
down electron tunnels into a state singly occupied by a 
spin-up electron, they form a pair; the interaction mixes 
this pair with the empty states at energies > Ez/2; the 
mixing is resonant at the energy E + and this resonance 
manifests itself as an anomalous contribution to the DoS. 
The position of the anomaly was predicted to be "uni- 
versal" , i.e. independent of both the dimensionality and 
the direction of the magnetic field. Recent experiments 
on superconducting Al films^ are in disagreement with 
the latter result: the position of the measured minimum 
in the DoS moves to higher energies as the component of 
the magnetic field perpendicular to the film is increased. 

In this paper we reconsider the calculation of the su- 
perconducting fluctuations corrections to the DoS in the 
normal phase for disordered films and wires in strong 
magnetic fields. To understand why the reconsideration 
is necessary, let us briefly review the qualitative argu- 
ment given in Ref. y to explain the position and width 
of the anomaly in films and wires: as discussed above, 
the position of the singularity was found to be located at 
E+ for tunneling into a grain; then a self-consistent ar- 
gument was given to find the characteristic energy scale 



Wd (width of the singularity in d dimensions) in one- and 
two-dimensional systems. The latter argument is based 
on the assumption that the system can be effectively di- 
vided into zero-dimensional patches whose size Lw d is 
then found self-consistently - this assumption however 
breaks down if dephasing happens on a scale shorter than 
the patches' size. In the presence of a perpendicular mag- 
netic field, this scale is given for films by the magnetic 
length Ih = (hc/eH) 1 / 2 , and the break-down happens at 



where L~w d is>2 



j w d 
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This suggests that, for strong enough magnetic fields, 
additional contributions to the DoS which are particu- 
lar to one- and two-dimensional systems might become 
relevant^ below we show that this is indeed the case and 
for films in a tilted field, we present the detailed deriva- 
tion of the dependence of the position of the singularity 
on the perpendicular component of the magnetic field. 
From now on, we use units with K = 1. 

To describe the superconducting systems, we consider 
the pairing Hamiltonian 7i = H.q + Hint , where the non- 
interacting part 7io is given by the sum of the kinetic 
energy, the Zeeman energy and the disorder potential, 
and Hi n t is, in second-quantized notation (g > is the 
coupling constant): 



n 



■nt = -g J d d r^\{r)^\{r)^ l {r)^{r) . 



The DoS is related to the imaginary part of the one parti- 
cle Green's function and the latter can be calculated with 
e.g. the diagrammatic technique^ As this procedure is 
standard, we skip intermediate steps and we simply quote 
the final answer— for the one-loop fluctuation correction 
to the DoS for spin down electrons at zero temperature 
in d-dimensions: 



5v d {e) 1 
= Im 

U 7T 
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where vq is bare DoS per spin, C is the Cooperon— 
C(e,q) = - ' 



■Dq 2 



(6) 
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with D the diffusion constant, and A is the fluctuation 
propagator: 



A(w,q) = — 



In 



E% + (-iu + Dq 2 )' 
A 5 



(7) 



As noticed in Ref. |jj due to its analytical properties this 
propagator can be separated into a "pole" part and a 
"cut" part: 



A = K p + A c 



where 



A p (u,,q) = 



2 A 2 



v 2Vt -i{uo - 0) + Dq 2 



(8) 
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and A c is implicitly defined by Eqs. ©-©■ We can ac- 
cordingly write 8v,i as a sum of two terms: <5^( e ) = 
(5^(e) 4- <5z/j(e). The contribution <5z/j(e) can be found 
in Ref. 4 and for example in d = 2 it is proportional to 
ln[ln(e — Ez)]', we give it here no further consideration, 
since this contribution is less divergent than the ones we 
calculate below. We note that the separation JSJ) is pos- 
sible provided that Ez > A; if this condition is satisfied, 
our results are applicable even for fields smaller than the 
parallel critical field if c y as long as the sample is in the 
normal state. 

According to the above argument, we want to evaluate 
the right hand side of Eq. (JSJ with the substitution 



A ^ A p 



(10) 



and the result of this calculation gives the lowest order 
perturbative correction to the DoS. On the other hand 
in Ref. a resummation of perturbation theory was per- 
formed in the long wavelength approximation Dq 2 < A, 
but this approximation is not always applicable. Indeed 
let us consider a 2D system, in which case the approxi- 
mation means that, in energy units, the exchanged mo- 
mentum is limited by A/<?, where g — Atti/qD is the adi- 
mensional conductance; for good conductors A/g A 
and this energy scale is much smaller than the gap A. 
Since the Cooper pair energy [i.e. the position of the 
pole in Eq. 0] £1 is > A and the exchanged energy 
is much smaller than f2, we can effectively neglect the 
"Fermi sea" under the pair: from a formal point of view, 
we can extend the upper limit of integration in Eq. (j3J) 
to infinity. In other words, in this approximation we can 
neglect the exclusion principle, which forbids the elec- 
trons from "going under the sea" , i.e. interacting with 
electrons at energies below the Fermi energy^S However 
when a magnetic field with a component perpendicular 
to the film is present, Landau quantization renders the 
momentum variable discrete: 



Dq 2 -^n H [n 



where 



n H = AeDH sin0 /c 



(11) 
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FIG. 1: Width and position of the minimum for small (left) 
and large (right) perpendicular field as explained in the text. 



is the Cooperon cyclotron frequency and 9 is the tilt- 
ing angle. The momentum transfer (in energy units) is 
therefore of order Qh and when this becomes sufficiently 
large we are not allowed to neglect the exclusion prin- 
ciple anymore: the correct limits of integration must be 
considered^ For strong fields, the effect of the "Fermi 
sea" is to move the position of the minimum to higher 
energies*^ for <C Wi the anomaly, centered at E + 
has a width W2 much smaller than E + [cf. Eqs. and 
(|18fl ] and the Fermi sea is not probed by the excitations 
contributing to the anomaly, see Fig. On the con- 
trary, for Qh ^> W2 a sizable fraction of these excitations 
would be located below the Fermi energy, as shown by 
the dashed line in Fig.^D. However the exclusion princi- 
ple suppresses the contributions from under the "Fermi 
sea" , the profile of the anomaly becomes asymmetric and 
as a consequence the minimum appears to shift to higher 
energies, as in the solid line of Fig. QJd. This qualita- 
tive argument predicts that as we increase f2#, the shift 
also increases, in agreement with the quantitative result 
derived below. It also places the transition between the 
weak and strong field regimes at £Ih ~ W2 , which is the 
same condition as in Eq. @ . In summary, the long wave- 
length approximation is justified in the weak field regime 
and the non-perturbative approach of Ref. is necessary 
in this case; in the strong field limit, on the contrary, 
a perturbative calculation is sufficient, as we explicitly 
show below. 

After the substitution l|10(l . the integration over uo in 
Eq. JSJ can be performed exactly and we write the result 
as: 



(13) 



where we separated different contributions based on their 
degree of divergence in the parallel field case: 8v™ gives 
the most divergent term (a (^-function in d = 2) which was 
considered in Ref. 0, along with subleading terms; other 
subleading terms are collected in the the (less) divergent 
part 5v l d ; finally 8vl contains only finite contributions 
and hence we do not need its explicit form. The relevant 
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contributions are: 
&vf{e) _ A : 
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C(-Q,g) 
^C(e-£+,g) (14b) 
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x[C(e-E z ,q)-C(2e-E z ,q)} 
-— ^- [C(2e - £ z , g) - C(-fi, g 



with £4. and defined in Eqs. and respectively 
and we introduced: 



= (E z - n) /2 . 



(15) 



Note that Eq. I|14a|l reduces to the perturbative result 
of Ref. 3 upon replacing the square bracket with —2ni. 
In what follows we retain only the most divergent terms, 
as finite and weakly divergent contributions have already 
been discarded in performing the replacement (JTUJ. 

We now restrict our attention to the case d = 2 with a 
non-zero perpendicular component of the magnetic field, 
so that we must use the substitution (jllfl and replace in 
Eqs. i|14fl the integrations over momentum with summa- 
tions: 



d 2 q 

W) 2 



AttD 
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(16) 



n=0 
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The result for the most divergent contribution is: 
/ 

+ / 
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where the energy 



W 2 = A 2 /4 g n 



(18) 



characterizes the width of the tunneling anomaly in the 
parallel field and the function f(a; b) is defined as: 



/(a; b) = Im 



" dt_ 



it — ia 



where ip' is the derivative of the digamma function. The 
validity of Eq. (|17fl is restricted to the strong field regime: 



0.000 




FIG. 2: Tunneling anomaly in tilted field. The continuous 
lines are the correction 8u^ of Eq. 1171 for 9 — 35°, 23° and 
16° (top to bottom) with g = 5, Ez = 2.5A and gL = 2. 
The vertical dotted line is at e — E+. The dashed line is the 
approximate contribution of Ref. for 9 = 35° , which can 
also be obtained by letting b — *■ oo in Eq. 11191 . 



minimum is shifted to a higher energy; the contribution 
at small energies (near the Fermi sea) is suppressed; the 
overall shape is asymmetric about the minimum, in qual- 
itative agreement with the experiments^^ In addition, 
the anomaly is smaller than the prediction of the approx- 
imate formula - this could be relevant in a quantitative 
comparison with experiments. 

For the purpose of calculating the field-dependent po- 
sition of the minimum, we can differentiate Eq. (|17|) and 
then expand the result to first order in 5 = (e — E + )/Qh 
by assuming 8 -C h. Performing this calculation we find: 



n 2 



(21) 



where 

(3 = - 



<2tt 2 

with the prime denoting the derivative with respect to 
the first argument and the dot with respect to the second 
argument of /, and 
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H 



(23) 



n H > 7r 2 w 2 



(20) 



in which the correction is indeed smaller than the bare 
DoS. As an example, we plot both Eq. l(T7|) and the cor- 
responding approximate formula of Ref. |3j in Fig. [21 We 
note that, in agreement with our previous discussion: the 



While the above inequality is exact (and the equality is 
valid for Q/D,jj = 0), the one in Eq. should be un- 
derstood as follows: for E-,Sl > f2_ff/3 the approximate 
equality holds, while at lower values of these parameters 
is the upper limit that applies. Setting Eq. (|21|l to and 
using for an estimate the limiting values of (3 and 7, we 
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find for the position E* of the minimum: 

E* = E + +an H , a=-^"Qj - 0.055, (24) 

which generalizes the result of Ref. and reduces to it 
for fifl- = 0»i£ The smallness of the constant a justifies 
a posteriori the expansion. In the experiment of Ref. |6J, 
Eq. 12411 has been successfully tested; below we comment 
on the applicability of this equation for experimentally 
relevant values of the parameters and we give more reli- 
able estimates for the dependence of E* on fl^. 

According to the definitions Q), @ and l|15|l . the in- 
equalities E- < il < E + < Ez hold; in the experiments 6 
all these quantities are ~ A. This means that Wi ~ 
A/4g -C A for good conductors. On the other hand 

9jL = ± gsin () > (25) 
Ez 9l 

so that for "large" tilting angles we have Vt h 3> Ez\ thus 
we find that the condition (|20|l is experimentally satis- 
fied. However the same reasoning shows that the con- 
ditions for the applicability of the upper limit estimate 
in Eq. (1221) are easily violated. A more detailed study 
of the correction l|17ll shows that for the position E* of 
the minimum we can distinguish two regimes for different 
ranges of the parameter £Ih/Ez, namely: 

E*=E ++ a,(^,^a H , l<g <7 (26a) 

E* = 0.543S+ + aoo n H , r-r- > 7 (26b) 

Ez 

with 

aoo~ 0.144. (27) 

The coefficient ot\ depends weakly on the field through 
Ez] for fields larger than about twice the parallel critical 
field this dependence can be neglected^ and we find: 

ax = ||^ (2) Q) W Q) * 0M6±gsin9 (28) 

independent of the field. The result l|26al) is obtained 
by expanding (3 [Eq. (|22|) ] as a function of the quantities 
E_/tt H , fl/Q H , while Eq. (|26bp is found by taking the 
limit g — > oo in the derivative of Eq. (|17|l : in this case the 



finite value of 5 [defined before Eq. J2J] is calculated nu- 
merically and then we evaluate the coefficients of the first 
order corrections in E-/£Ih, VI/Q,h - these are respon- 
sible for the suppression of the E + term. We stress that 
in the tilted field the linear dependence of the position 
of the minimum on the field is a robust prediction. We 
also notice that in the allowed region of the parameters, 
the coefficient ct\ in front of Qh varies between 0.016 and 
. 1 1 2 , i . e . it agrees with a of Eq. I|24|) within a factor of 3 ; 
since it is also true that a^/a < 3, the estimate Eq. l(24*|) 
is a correct order-of-magnitude one. Indeed that equation 
has been successfully applied in a study of a g = 5.6 sam- 
ple for whichi^ VLh/Ez < 8 in Ref. Experiments in 
higher conductance samples, where D,h/Ez ^> 7, would 
enable to test Eq. I|26b|) directly: its simpler dependence 
on the physical parameters as compared to Eq. I|26a|) 
would make possible a more quantitative comparison be- 
tween theory and experiments. 

Up to now, we have not considered the phase relaxation 
effect of the (parallel component of the) magnetic field, 
since it can be usually neglected in the tilted fields This 
effect is accounted for by shifting the frequency in the 
definition of the Cooperon: e — > e + i/th, with* 

1/th = n 2 H{l /48E T . (29) 

Here f2#n is as in Eq. (|T2"|) but with cos 9 instead of 
sin 9 and the transverse Thouless energy is: Et = D/a 2 , 
where a is the film's width. In the parallel field, mo- 
mentum integration in Eqs. (|14|) is straightforward; we 
do not give here the explicit answers (due to space lim- 
itations), but the qualitative results are the same as for 
the tilted field upon replacing £Ih —> ^/t~h- Also, we 
note that Eqs. (|14|) can be used to evaluate the correc- 
tion to the DoS of superconducting wires; for both wires 
and the parallel field case, the one-loop approximation is 
valid when Warn 

In conclusion, we calculated the one-loop interaction 
correction to the Density of States for superconducting 
films in tilted field and in the paramagnetic phase. We 
found a minimum in the Density of States which shifts 
linearly with the applied field to higher energies, see 
e.g. Eq. I|24|) : this dependence has been demonstrated 
experimentally. 6 
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